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Abstract 

We present an ab initio study of the interface energies of cubic yttria-stabilized zirconia (YSZ) 
epitaxial layers on a (0001)q,_ai2O3 substrate. The interfaces are modelled using a supercell ge- 
ometry and the calculations are carried out in the framework of density-functional theory (DFT) 
and the local-density approximation (LDA) using the projector-augmented-wave (PAW) pseudopo- 
tential approach. Our calculations clearly demonstrate that the (111)ysz||(0001)q,_ai2O3 interface 
energy is lower than that of (100)ysz||(0001)q,_ai2O3- This result is central to understanding the 
behaviour of YSZ thin solid film islanding on (0001)ci,_ai2O3 substrates, either flat or in presence 
of defects. 

PACS numbers: 68.35.bt, 68.55.aj 
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I. INTRODUCTION 



Over the last few years, many experimental efforts have been expended on the fabrication 
of self-patterned, epitaxial nanocrystals (metallic, semiconductor or oxyde) on crystalline 
substrates.-'^ii^i^'^'^ The aim is to synthesize homogeneous patterns of epitaxial crystals in 
order to induce quantum confinement — intimately related to the shape and the size of the 
nanocrystals — in order to achieve enhanced optical and/or magnetic properties.->^'2i'i2i»ii'i2i'i^ 
Several theoretical investigations have been concerned with the physical parameters respon- 
sible for the geometric properties of the nanocrystals.— 

The fabrication of such systems can be realised through various techniques involving the 
formation of nanometer-scale islands in a collective way, a process known as self-organisation. 
The basic idea is to promote the formation of nanocrystals on a crystalline substrate using 
thin solid films which demonstrate spontaneous evolution from a continuous 2D solid layer 
to a rough and/or discontinuous film, i.e., 3D epitaxial nanocrystals. The formation of the 
nanocrystals takes place during or after the deposition of the film as a way of reducing 
the total energy of the epitaxial {layer 1 1 substrate} system by the relaxation of the inter- 
face and/or surface stresses and strains. We briefly describe, in what follows, three of the 
most popular "bottom-up" approaches involved in epitaxial nanocrystals self-organisation 
processes. 

A first approach is chemical vapor deposition (CVD), which consists in mixing chemical 
species in vacuum, which then react or decompose on the surface of a substrate to form 
a thin solid film. A second approach is physical vapor deposition (PVD), whereby matter 
is extracted from a solid target with, for example, a laser or an ion beam; the extracted 
ions attach to the surface of a substrate and eventually constitute a thin solid film. CVD 
and PVD, which lead to the formation of nanocrystals during the deposition process, have 
been successfully applied to the fabrication of self-organized arrays of semiconductor or 
metal nanocrystals, commonly called quantum dots (QDs). One of the most widely studied 
system is {Si||Ge}; this is characterized by the formation of epitaxial, faceted Ge QDs either 
at the top of a continuous Ge wetting layer (Stranski-Krastanov growth), or directly at 
the surface of the Si substrate (Volmer- Weber growth) .—i^i2Si2^i2ii2^i2Si2I A third approach 
goes by the deposition of a continuous thin xerogel film at the surface of a substrate by 
sol-gel dip-coating.— Through thermal treatment, the continuous thin solid film crystallizes 
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and breaks into several crystals through surface diffusion,— i^SiiSi leading to the formation 
of discrete epitaxial islands on the surface of the substrate. In this case, unlike PVD or 
CVD, the formation of the epitaxial islands takes place after the deposition of the film. This 
technology is particularly efficient for designing arrays of oxide nano-islands, as recently 
demonstrated by Bachelet.— 

In this article, we are concerned with the {YSZ| |(0001)q,_ai2O3} system, where epitaxial 
nano-islands of cubic yttria-stabilized zirconia (YSZ) form during the thermal treatment of 
a thin xerogel film deposited by sol-gel dip-coating on a (0001)q,_ai2O3 substrate. Exper- 
imentally, Bachelet et al. have demonstrated that the shape and size of the YSZ islands 
are directly linked to their epitaxial relation with the substrate.— Indeed, for a substrate 
without defects, the islands are top-fiat with large interface areas and exhibit the following 
in-plane and out-of-plane crystallographic orientations: 

(100)ysz||(0001),_ai2O3, [001]ysz||[010]„_ai2O3, (1) 
(100)Ysz||(0001)„_AbO3, [001]ysz||[110]„-ai2O3, (2) 

whereas for a substrate containing defects, some islands are round and thicker than the 
top-fiat ones but with lower interface areas, and possess the in-plane and-out-of plane ori- 
entations: 

(111)ysz||(0001)«_ai2O3, [110]Ysz||[110]„-AbO3- (3) 

The morphological evolution of the thin solid film into discrete nano-islands proceeds by 
an abnormal grain growth driven by the interface during thermal treatment.— From our 
previous theoretical investigations of this system,'^ and in good agreement with experimental 
results, we demonstrated that the shape and size transition from top-fiat to round is linked 
to the presence of defects at the surface of the substrate which induce enhanced growth in 
height. 

However, the preferred formation of interfaces (1) and (2) over interface (3) for a perfect 
(0001)a_Ai2O3 substrate is not clearly understood. On the basis of energy considerations it 
can be argued that, to first order, an epitaxial crystal is in equilibrium with both the vacuum 
through the free surface energy and with the substrate through the interface energy, where 
both energies are related to the crystallographic orientations. Using ab initio methods, 
Ballabio et al. have demonstrated that the free surface energy of (IOO)ysz is higher than 
that of (III)ysz-— 
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In this article, we propose to examine the interface energies defined by the epitaxial 
relations (1), (2), and (3) on a perfect (0001)a_Ai2O3 substrate. We argue that the knowl- 
edge of the interface energies is sufficient for a proper comparison of the behaviour of the 
three interfaces. There are several theoretical investigations of {metal| |oxide} interfaces, 
in particular for (0001)q,_ai2O3 because of its technological significance in thermal barrier 
coatings and catalytic devices.— '^'^'^'^i^'^'^ However, few studies have been devoted to 
{oxide| |oxide} interfaces^ and, to the best of our knowledge, none have been concerned with 
the {YSZ||(0001)„ -AI2O3} system. Our calculations clearly demonstrate that interface (3) 
is energetically favored over interfaces (1) and (2). We therefore propose a general expla- 
nation for the behaviour of the islanding process of YSZ thin solid films on (0001)q,_ai2O3 
substrates, either perfect or with surface defects, in the light of experimental and theoretical 
invest igat ions .—^ 

The article is constructed as follows. In section [TTl we give the details of the numerical 
procedure, followed in section IIIII by a demonstration of the ability of the PAW pseudopo- 
tentials to reproduce the correct structural properties of a and fi;-Al203, the low-pressure 
polymorphs of Zr02, and the Y2O3 bixbyite structure. In section [IV] we first discuss the 
calculation of the unrelaxed and relaxed stoichiometric free surface energies, then we present 
the atomic-scale models for the the (1), (2) and (3) {YSZ||(0001)a_Ai2O3} interfaces and the 
results of our calculations. A general conclusion is provided in Section |Vl 

II. COMPUTATIONAL DETAILS 

All calculations were carried out in the framework of density functional theory (DFT) 
using the Abinit code,^"^ where the wave functions are expanded in plane waves. The atomic 
pseudopotentials were constructed with the atompaw program^^ within the frozen-core ap- 
proximation, using the projector-augmented-waves (PAW) method originally proposed by 
Blochl.^^ For the exchange- correlation functional, we employed the local density approxi- 
mation (LDA) as parametrized by Perdew and Wang.— The atomic wave functions were 
augmented with 3, 6, 5, and 5 projectors within a spherical augmentation region of radii 
1.4, 1.8, 2, and 2 Bohrs for O, Al, Y, and Zr atoms, respectively. The 2s and 3s semi-core 
states of Al, as well as the 4s and 4p semi-core states of Y and Zr, were treated as va- 
lence states to generate the pseudo-wave and projector functions within the augmentation 
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region. We found that taking Y and Zr semi-core states as valence is debatable. Indeed, 
Jansen^ demonstrated that, due to the large energy difference between O and Zr(4s,4p) 
resonances, the 0-Zr(4s,4j9) hybridization is weak in Zr02. This argument was applied 
by Christensen and Carter to study the free surfaces of Zr02 low-pressure polymorphs^ 
and the {{001)zrO2\\{^^'^'^)a-Ai2O3} interface.'^^ However, in previous studies of bulk YSZ by 
Stapper et aL^^ and of YSZ slabs by Ballabio et al.,— the (4s, 4p) semi-core states of Y and 
Zr were treated as valence. Here, we are also dealing with YSZ bulk and slab structures, as 
it is our purpose to characterize the {YSZ||(0001)q_ai2O3} interfaces. In section UTT] we will 
demonstrate that the treatment of 4s and Ap semi-core states of Y and Zr as valence states 
is appropriate to accurately describe the structural properties of Y2O3 bixbyite and Zr02 
low-pressure polymorphs. 

Additional details are as follows. For the a and n phases of AI2O3, as well as the Zr02 
low-pressure polymorphs and the Y2O3 bixbyite structure, all discussed in section llllt we 
used a kinetic energy cutoff of 15 Ha (~ 408 eV) and 2x2x2 Monkhorst-Pack grid^ 
for the Brillouin-zone integrations of bulk unit cells, which is a standard choice for wide 
band gap oxides.— With these parameters, the total energies are converged to within 10~^ 
Ha/atom (0.2 eV/atom) and the forces to better than 10~^ Ha/ (Bohr. atom). The atoms 
were relaxed to their ground-state positions using the Broyden-Fletcher-Goldfarb-Shanno 
(BFGS) algorithm..^ 

The slab geometry for the (0001)a_Ai2O3 system is described in section HVl in this case 
we used a 2 x 2 x 1 Monkhorst-Pack grid in the z direction. For the (1), (2) and (3) 
{YSZ||(0001)q,_ai2O3} interfaces, and according to previous studies,—!^ only the P point 
was used to integrate the Brillouin zone. The convergence criteria for relaxation were the 
same as above. 

III. STRUCTURAL PARAMETERS OF THE BULK PHASES 

We present here our results for the various bulk phases in order to ensure that our ap- 
proach yields the correct structural parameters. For this purpose, we computed the relaxed 
lattice parameters and ionic positions (we follow the Wyckoff convention^'^) for each crys- 
talline structure. We also computed the ground state energies E as a. function of the volume 
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FIG. 1: (Color online) Energy vs volume for a-Al203 and K-AI2O3, as indicated 
V of the unit cell; these can be fitted to the Murnaghan equation of state:— 

with 

thus yielding the equilibrium volume Vq and the bulk modulus P is the pressure. 

Q;-Al203 (also known as corundum) is the thermodynamically stable polymorph of alu- 
mina at low pressure. It is widely used for epitaxial thin solid film growth, in particular the 
(0001) "C-cut" and (1012) "R-cut" families of planes. It has space group i?3c and can be 
represented either by a rhombohedral unit cell with 10 atoms (2 AI2O3 molecular units) or a 
hexagonal cell with 30 atoms (6 AI2O3 molecular units). In our calculations, we assigned the 
initial positions of the Al and O ions according to the (experimental) values of Wyckoff.— In 
the conventional hexagonal cell, there are 6 oxygen planes organised in the ...ABABAB... 
closed-packed stacking sequence in the c direction, with the aluminum ions occupying 2/3 
of the octahedral sites. 

Apart from the a phase, there are two metastable polymorphs of alumina which are of 
practical interest, viz. K-AI2O3 and 7 — AI2O3; we focus on K-AI2O3 hereafter. This poly- 
morph can be synthesized by CVD and, because of its hardness, is used as surface coating.— 
The crystalline structure of K-AI2O3 is orthorhombic, space group Pna2i, with 40 atoms (8 
AI2O3 molecular units) in the unit cell.— There are four oxygen planes in the cell, organized 
in the ...ABC ABC ABC... stacking sequence, and the aluminum ions fill the tetrahedral 
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TABLE I: Parameters of the Murnaghan equation of state, Vq (A'^/atom), Bq (GPa), and Bq, for 
a-Al203 and K-AI2O3, and comparison with other results from the hterature. 
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'Emp.' stands for 'empirical model' 



and/or octahedral sites.— The positions of the Al atoms in this structure is under debate. 
Indeed, Belonoshko et al.^^ proposed a model in which 2/3 of the octahedral sites are filled 
with Al ions as in a-Al203 whereas, according to the theoretical studies of Yourdshahyan 
et al., the most stable structure is one of the nine possible configurations for which the 
Al are only in octahedral positions.—"^ In this work, for the sake of simplicity and clar- 
ity, we have chosen to fix the reduced coordinates of Al and O ions to the experimental 
values issued from the Rietveld refinement of Smrcok and al.;^^ thus, all the ions are in po- 
sition 4a with coordinates {xi,yi, zi) for Al and (x2, 2/2, -22) for O, with symmetry operations 
(x, y, z- -X, -y, z + 1/2; x + 1/2, -y + 1/2, z- -x + 1/2, y + 1/2, z + 1/2). 

The Eiy") curves for both polymorphs are presented in Fig. [T] and the parameters of 
the Murnaghan equation of state are provided in Table [B One can see that our ah initio 
calculations do reproduce the correct relative stability of the two phases; it is indeed known 
that at low pressure the fi;-Al203 — > Q;-Al203 phase transition occurs around 1000°C. 
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FIG. 2: (Color online) Energy vs volume for the low-pressure polymorphs of Zr02, as indicated 

The lattice parameters of the rhombohedral and hexagonal phases of a-Al203 are sum- 
marized in Table HTl The parameters for the hexagonal cell, a^i, 0,^2, and Ch, are deduced 
from the rhombohedral ones, cirhi, arh2, and cirhs, as follows:— 

(^hl = (^vhl — (^rh2 (6) 
3.h2 = arh2 — flrhS (7) 
Ch = flrhl + 5rh2 + ^rhS (8) 

with Orh = llflrhill = ||arh2|| = llflrhsll, o,h = \\(^hi\\ = ||f?h2||, and Ch = \\ch\\- The parameters 
for (orthorhombic) fi;-Al203 are Oq, &o and Cq. 

We now turn to zirconium dioxide, an important material for applications in optical, 
mechanical and thermal coatings. Here we are concerned with the low-pressure polymorphs 
of Zr02 stoichiometry. From to 1400 K, Zr02 is monoclinic ('m'; this phase is called 
baddeleyite) , of space group P2i/cf^ between 1400 and 2650 K, it is tetragonal ('t'), of 
space group P42/nmc;— finally, above 2650 K and all the way to the melting point, it 
is cubic ('c'), of space group Fm3m.— Here we describe the three phases through their 
conventional unit cells with 12 atoms (4 Zr02 molecular units). Table IIIII presents the 
parameters of the Murnaghan equation of state derived from the E{y) curves of Fig. [2l 
Our calculations reproduce the correct relative stability of the three polymorphs and our 
fitted parameters agree with previous theoretical and experimental results. The structural 
parameters are provided in Table IIVI 

Finally, we discuss the structural parameters of Y2O3 (bixbyite), which is body-centered 
cubic, space group Ia3>. This can be viewed in the conventional unit cell with 80 atoms 
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TABLE II: Lattice parameters for a-Al203 in the rhombohedral and hexagonal cehs, and for k- 
AI2O3 in the orthorhombic cell (A), and comparison with other results from the literature; is 
the angle of the rhombohedral cell (degrees). 
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and lattice parameter a = 10.604 A,^^ or in the primitive unit cell with 40 atoms and 
lattice parameter a' = a-\/3/2. In this work, we have used the primitive unit cell; there 
are two nonequivalent yttrium sites, 8a for Yj of coordinates (1/4,1/4,1/4) and 24(1 for 
Y// of coordinates (m,0, 1/4). The oxygen site is 48e of coordinates {x,y,z). The relevant 
symmetry operations can be found elsewhere.— Table |V] presents the parameters of the 
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TABLE III: Parameters of the Murnag 
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Murnaghan equation of state, the lattice parameter, and the internal ionic positions for Y// 
and O. 

The above results clearly establish the ability of the PAW method to reproduce the correct 
structural properties of the systems we are concerned with. In what follows, we present first 
the models used to simulate the (0001) surface of a-Al203 as well as the (100) and (111) 
surfaces of YSZ. We then discuss the method for constructing the interface supercells, which 
will be used to calculate the interface energies. 
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TABLE IV: Lattice parameters for m-, t-, and c-Zr02 (A), and comparison with other results from 
the hterature; (3 is the angle of the monoclinic phase (degrees) and dz is the tetragonal distorsion 
of O atoms in the c direction of the tetragonal phase. 
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IV. INTERFACE ENERGIES OF {YSZ||(0001)c,-Ai2O3} (1), (2) AND (3) 

In order to calculate the interface energies, we constructed supercell models consisting 
of a (0001)q,_ai2O3 slab for the substrate and a (IOO)ysz or (III)ysz slab for the epitaxial 
layer. The supercells are parallelepipeds, and the interface is taken to be perpendicular to 
the z direction. Before proceeding, however, we consider the free surfaces and compute their 
unrelaxed and relaxed energies. Periodic boundary conditions are used; for free surfaces, a 
vacuum region is inserted in the supercell. Thus, in all cases there are two interfaces, either 
between the two materials or between the surface of the material and the vacuum. 

The thickness of the slab must be sufficient to yield converged results and yet remain 
computationally manageable. To this end, one may first define two structures — one for the 
bulk and one for the slab — having the same number of atomic layers A^; the slab has two 
free surfaces owing to the presence of a vacuum region (see above). The free surface energy 
is then defined as the excess energy of the slab [{hkl) indices] relative to the bulk divided 
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TABLE V: Lattice parameter a (A), parameters of the Murnaghan equation of state, Vq (A'^/atom), 
Bq (GPa) and B'q, and Wyckoff coordinates of Yjj and O for the Y2O3 bixbyite structure, and 
comparison with other results from the hterature. 
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TABLE VI: Free surface energies 7(oooi) (J/m^) of the unrelaxed and relaxed a-Al203 models for 
iV = 9, 12, 15, 18 atomic layers (i.e., AN = 3 here). 
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9 
12 
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^slab(iV)-^AEbulk(iV) 

s 

AEbuik(iV) = ^buik(iV) - ^buik(iV - AN), (10) 



IhkiiN) = (9) 



where AA^ is the difference in the number of layers between two different slab models. The 
convergence of 'jhki can be studied as a function of A^. 

However, Boettger has shown that this approach is not very accurate,— as the convergence 

12 



of 'jhki{N) depends on the thickness of both the slab and the bulk, but also on Ai?buik(^)- 
He proposed to use, instead: 

lhki[J^) = (llj 

AE,i,b(iV') = ^siab(iV') - ^siab(iV' - AN) (12) 

where N' is the number of layers for which the value of Aii^siab is sufficiently converged to 
ensure that the behaviour of '^hki is a function of only. Thus, in this approach, the term 
Ai?buik(A^) in Eq- Elis replaced by Ai?siab(A^') in Eq. [HI thereby reducing the convergence 
study of to slab calculations. For the sake of comparison with previous ab initio results 
on this system,— we used the Boettger method here. 



A. (0001)q_ai2O3 free surface energy 

Several calculations of stoichiometric and non stoichiometric free surface energies of 
(0001)a-Ai2O3 have been reported in the literature.— '^i^'^'^'^'^ Here we deal only with 
stoichiometric systems; this choice is not restrictive as it was demonstrated that the most 
stable (0001) surface of a-Al203 is stoichiometric, Al-terminated, in a wide range of Pq^.—^ 
Further, previous ab initio calculations^ have shown that a vacuum thickness of 10 A is suf- 
ficient and this is the value we have used. We have nevertheless studied the convergence 
with regard to the thickness of the solid, viz. 9, 12, 15, and 18 atomic layers. 

The results, presented in Table IVIl are found to be in very good agreement with those 
of Siegel.— One may note the huge differences between the unrelaxed and relaxed energies 
— the absolute differences are ~2 J/m^. This is a consequence of the inward relaxation of 
the atomic planes in the z direction. Table IVIII gives the average relaxation of the atomic 
planes relative to the original bulk spacing. As found in previous calculations, the inward 
relaxation of the Al atomic plane is close to 80% and leads to the formation of sp^-like 
atomic bonding at the free surface. The inward relaxation of the top Al plane is related 
to the increase of the electronic density, yielding a lower free surface energy. Our study 
demonstrates that 15 atomic layers are needed to model accurately the bulk structure, but 
the results are already quite reasonable for = 9, offering a good compromise between 
accuracy and computational workload as we discuss in Sec. IIVCI 
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TABLE VII: Average relaxation of the atomic planes in the [001] direction for a-Al2 03, expressed 
as a proportion of the initial bulk spacing for = 9, 12, 15, 18 atomic layers. 
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B. (IOO)ysz and (III)ysz free surface energies 

We now consider the YSZ (100) and (111) free surface energies. We follow the approach 
proposed by Stapper et al— and Ballabio et al.^^ to build the bulk and slab structures. YSZ 
is a solid solution of Y2O3 in Zr02, of space group FmSm (as c-Zr02). Proper simulation 
of YSZ depends on two parameters: (i) The size of the simulation cell, which must be large 
enough to provide a good statistical representation of the proportion of Y atoms and O 
vacancies (Vq) for a given molar proportion of Y2O3. (ii) The ground state energy, which 
depends on the relative positions of Y ions and O vacancies, and which cannot be chosen 
at random: Stapper et al. have indeed shown that the most stable configuration is that for 
which the O vacancies are next-nearest neighbours to yttrium atoms.— 

Here, the doping level of Y2O3 is set to 10% molar, consistent with the experimental 
studies of Bachelet et alM For consistency and comparison with previous works, the positions 
of the Y ions and the O vacancies are chosen such that two Vq's cannot be closer to one 
another than third nearest neighbour; two Y's can be nearest neighbours but a Y cannot 
be closer to a Vb than next nearest neighbour. The bulk cell of (IOO)ysz is made up of 4 
(Zr,Y) and 5 (0,Vb) atomic layers in the [100] direction (A^ = 9), for a total of 93 atoms 
(26 Zr, 61 O, 6 Y; 3 Vq). The dimensions of the bulk cell are 6x = 6y = 6z = 2aYsz, where 
the theoretical lattice parameter is derived from the experimental relation established by 



14 



TABLE VIII: Unrelaxed and relaxed stoichiometric free surface energies 7(ioo) and 7{iii) for YSZ 
(J/m2). 



Unrelaxed 



Relaxed 



7(100) 
2.79 



7(111) 
1.30 " 



7(100) 
1.71 

1.75 



T^(iii) 
1.17 " 



1.04 



'^Thisrwprk 
^Ref. y 



Pascual and Duran,*^ aysz = c^o + 0.003x, with x the molar percent of Y2O3. Using the 
values of oq given in Table HVknd setting x = 0.1, we obtain aysz = 5.05 A. For the (III)ysz 
cell, we have 3 (Zr,Y) and 6 (0,Vo) atomic layers in the [111] direction (A^ = 9), for a 
total of 140 atoms (40 Zr, 92 O, 8 Y; 4 Vq). The dimensions of the cell are 5x = 2v^aYsz, 



Sy = 2A/372aYsz, and SzVSaysz- 

From these bulk cells, surface slabs are constructed by introducing a 10 A-thick vacuum 
layer along z. For consistency with the case of (Al-terminated) a-Al203, we also consider 
stoichiometric surfaces hereafter so that, in both cases, the cells are terminated by an oxygen 
plane. In the case of (IOO)ysz, this requires half of the O atoms to be removed from each side 
of the slab. We placed one oxygen vacancy on each side of the slab cells on the free surfaces. 
The (100) and (111) free surface energies of YSZ, computed using Eq. [9|, are presented in 
Table IVnH the areas of the free surfaces are S = Saygz and S = Sv^flysz; respectively. 

We find good agreement with Ballabio et al. for the relaxed value of 7(ioo)) but there is 
a small difference for 7(111); this might be due to the use of different relaxation schemes — 
BFGS here, vs Car-Parinello molecular dynamics®^ in Ref. 3^. More important, however, 
we observe large changes arising from relaxation: ~1.2 J/m^ for 7(100) and ~0.2 J/m^ for 
7(111), with an average inward relaxation of the top O plane of ~25% and ~8%, respectively. 
These results are not surprising since (111) corresponds to a dense arrangement of the atomic 
planes, which is not the case for (100). 
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FIG. 3: Relative crystallographic orientations and dimensions of the a-Al2 03 and YSZ unit cells 
(dashed lines) in the x and y directions for each interface model. For model (3), the YSZ [100] and 
[010] crystallographic orientations are out of the (x, y) plane and are represented as dotted lines. 



C. {YSZ||(0001)a-Ai2O3} interfaces 

We now turn to the (1), (2), and (3) {YSZ| |(0001)a_Ai2O3} interface models, constructed 
from the structures discussed in the previous sections. More specifically, several unit cells 
must be assembled in the (x,?/) plane so as to minimize the lattice mismatch; since the 
computational effort increases very rapidly with size, the number of cells of each material 
must be chosen such that the mismatch is no larger than a few percent in each supercell. 

Figure [3] provides a schematic representation of the relative orientations of the two ma- 
terials; ball-and-stick models are presented in Fig. HI In practice, keeping the mismatch to 
within a few percent would require unit cells containing at least 500 atoms. This is clearly 
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TABLE IX: Lattice mismatches e relative to the initial bulk spacings (%) and total number of 
atoms ria for each interface model. Ni^ and Ngx^ and Niy and Ngy^ are the number of unit cells for 
the layer and the substrate in the x and y directions, respectively. 



iV/x = 3 Niy = l 

= 2 

= 519 



Nsx =4 Nsy = 2 



eix = 3.93 eiy = -3.30 

(-SX = -3.64 tsv = 3.53 



(2) 



Nix = 1 
Nsx = 2 
€ix = -3.30 
^sx — 3.53 



(3) 



Nix = 1 
Nsx = 3 
eix = -0.45 
esx = 0.46 



Nly 



3 

: 4 



e^j^ = 3.93 



519 



-3.64 



Nly 



2 
3 



"sy - 

eiy = -0.46 



550 



0.47 



unmanageable. In order to reduce the workload, we fixed the thickness of the a-Al203 sub- 
strate to 9 atomic layers; as mentioned earlier, while the system parameters are not fully 
converged at this value, they are nevertheless adequately described. In addition, as we will 
be comparing different structures with the same number of layers, systematic errors will 
cancel out to a large extent. The computational workload remains considerable, but was 
alleviated by running the calculations in parallel on up to 252 processors. 

The details of the geometry of the three systems are given in Table llXi the mismatches 
are in all cases less than ~4%. The distance between the substrate and the layer was fixed 
by assuming that the YSZ O plane at the interface lies at the position where an a-Al203 O 
plane would have been in an infinite system. 

The total energy Ei^t of a given model structure can be related to the ideal work of 
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FIG. 4: (Color online) Side views of the unrelaxed interface models (1), (2) and (3). The Y and 
Zr atoms are white, the O atoms are light bluel&id the Al atoms are purple. The crystallographic 
orientations are those of the YSZ phase. 



(b) 
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FIG. 5: (Color online) Schematic illustration of the geometry of the system and corresponding in- 
terface energies: periodically-replicated model interface structure (a) with vacuum and (b) without 



vacuum. 



adhesion VTad as follows: 



ad 



Es + El — 
S 



(13) 



where Es and Ei are the total energies of the substrate and the layer, respectively, and S is 
the area of the interface. PVad can be expressed in terms of the interface and surface energies 
as: 

Vrad = 27, + 27/-7int-7int', (14) 

where 7^ = 7a-Ai203 = 7{oooi) and 7; = 7ysz = 7{ioo) for interfaces (1) and (2) = 7(111) for 
interface (3). The quantity 7int' is the energy resulting from the presence of two interfaces 
in absence of vacuum. The significance of the various quantities entering Eqs. [13] and [H] in 
relation to the geometry of the models is schematically illustrated in Fig. [51 

One may argue that the interface energy between two solids is lower than the sum of the 
free energies of the two surfaces: 



Combining Eqs. 



7int' < 7(0001) + 7ysz- 

m and [15] yields a lower bound to the interface energy: 



7int > 



Ejnt — iEs + El 

s 



7(0001) + 7ysz- 



(15) 



(16) 



This is useful for comparing interface energies as 7ysz and 7(oooi) are known. However, 
because the mismatches are finite (cf. Table [TX]) . the substrate and layers are both under 
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strain and, as a consequence, the values of Eg, Ei, 7ysz, and 7(oooi) must reflect this. In 
practice, the values are adjusted for the actual lattice parameters of the substrate and the 
layers; the corresponding strained, unrelaxed free surface energies are listed in Table |Xl 

One may note that for both layers and substrates, the strained unrelaxed free surface 
energies are strictly lower than their corresponding unstrained values (see Tables IVII and 
IVIIII for comparison). This behaviour is related to both the evolution of the free surface 
(contraction or dilatation, see Table HXl) and the modification of the electronic density. Using 
Eqs. [15] and [161 and the parameters listed in Table [Xj we have 

4.48 < 7i,t(i) < 6.81, (17) 
4.41 < 7int(2) < 6.89, (18) 
2.46 < 7i,t(3) < 5.36. (19) 

Our calculations demonstrate, therefore, that the interface energy for model (3) tends to 
be lower than that of models (1) and (2), implying that the fomer is the most thermody- 
namically stable of the three. This result is consistent with the continuity of the three-fold 
a — AI2O3 (0001) symmetry axis with the YSZ (111) symmetry axis in model (3),~ i.e., 
electronic bonding is much stronger here than in models (1) and (2). As a consequence, 
on a perfect a — AI2O3 substrate, the creation of a type-(3) interface requires much more 
energy than interfaces of types (1) and (2) because it is more strongly bonded. Thus, at the 
beginning of the islanding process, when the thin solid film still exhibits a large interface 
area with the substrate, the formation of interfaces (1) and (2) is favored over interface (3). 
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As a consequence, (III)ysz islands are not expected to form on a perfect a — AI2O3 substrate 
— only (IOO)ysz islands with large interface areas should be observed. 

On an imperfect substrate, now, the islands nucleating at the location of the defects 
are subject to enhanced growth in height.™ As a result, the interface area decreases and 
the energy cost required for the island to create interface (3) decreases with regard to the 
whole internal energy of the island. This allows one to understand why, on rough substrates, 
(III)ysz rounded islands are observed only at the location of defects. 

V. CONCLUSION 

We have demonstrated using an ab initio approach that the (111)ysz| |(0001)a_Ai2O3 inter- 
face is thermodynamically more stable than (100)ysz||(0001)q,_ai2O3 • This result allows us 
to understand and coherently describe the islanding process during the thermal treatment 
of YSZ on a — AI2O3, either perfect or with defects. On a perfect substrate, the formation 
of the (IOO)yszI |(0001)a_Ai2O3 interface is energetically favored over (111)ysz||(0001)q_ai2O35 
opening up the way to the formation of (IOO)ysz islands on the perfect substrate. On an im- 
perfect surface, the formation of islands at the location of defects leads to enhanced growth 
in height. As the interface area decreases, the energy cost required to form interface (3) 
does too. As a consequence, the (111)ysz||(0001)q_ai2O3 interface can form at the location 
of defects, which explains that both (IOO)ysz and (III)ysz islands are observed in this case. 
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